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Abstract 



> , 

^SJ , We consider the uncertainty of theoretical calculations for a specific difference of 

CN ■ the hyperfine intervals in the Is and 2s states in a light hydrogen-like atom. For 

a number of crucial radiative corrections the result for hydrogen atom and helium 



I ion appears as an extrapolation of the numerical data from medium to low Z. An 

' approach to a plausible estimation of the uncertainty is suggested using the example 

■5^ ■ of the difference D21 = 8£'hfs(2s) — £'hfs(ls)- 

Q^l Key words: Tests of quantum electrodynamics, Hyperfine structure 
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Studies of the hyperfine structure in light hydrogen-like atoms are of inter- 
est because of a possibility for a precision test of the bound state Quantum 
Electrodynamics (QED). However, in spite of a record accuracy achieved in 
experiments on the hyperfine interval in the ground state of hydrogen, deu- 
terium, tritium and helium-3 ion (see, e.g., [1]), we should acknowledge that 
such a test cannot be really successful until effects of the nuclear structure are 
known accurately enough. 

In fact, the bound state QED effects are lower than those of the nuclear 
structure. Obviously, the latter cannot be known well enough. However, there 
is a hopeful opportunity to perform such a test if we consider data of two 
measurements: for hyperfine intervals in the ground and metastable states 
(i.e. for the Is and the 2s states). 
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Combining the difference 



(1) 



we take advantage of a substantial cancellation of various contributions caused 
by short distance effects, since 



where ipnsi^) is a value of the Schrodinger- Coulomb wave function at origin 
for the ns state. 

The nuclear structure is such a short-distance effect and its leading contri- 
bution is proportional to |'0(O)|^. That is the crucial feature of the difference 
D21 that a complete cancellation of the leading nuclear structure term takes 
place. Meanwhile the higher order nuclear effects are suppressed by the factor 
of {ZaY and are well under control (see [2,3] for detail). 

Theoretical contributions are conventionally presented in terms of the so-called 
Fermi energy, which determines the hyperfine interval for the ground state. 
For the nuclear spin 1/2 it is of the form 

EF^^Z'a'm-^, (3) 
6 Hb 



which involves the nuclear magnetic moment and the Bohr magneton ij,b- 

Throughout the paper we apply units in which h = c = 1; m is the electron 
mass; M is the nuclear mass; Z is the nuclear change; and while all equations 
are presented for the energy (A£'), the numerical results are for the related 
frequency (AE/h); all fractional values are in units of Ep. 

The leading contribution to the difference D21 is 

D^^^^liZafEp (4) 



and the known corrections include higher order QED terms (of the third and 
fourth order in the expansion over a, Za and m/M) and higher order nuclear 
effects. We mainly follow [3] after correcting a misprint in the expression for the 
nuclear magnetic moment of the helion (the nucleus of the helium-3 ion). The 
correction slightly shifts the result of [3] (D*f °'^(3He+) = -1 190.067(63) kHz) 
to L)*'j*=°''(3He+) = -1 190.083(63) kHz (see,e.g., [4]). Other corrections, but 
one, have a marginal effect and will be considered elsewhere. 
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Here, we consider the most important theoretical issue which affects our early 
predictions for Z = 1,2. Before discussing it let us remind that a substan- 
tial part of the theoretical uncertainty in calculations of [3] came from the 
one-loop self energy contribution estimated according to Ref. [5]. It is of ap- 
proximately the same value as contributions to the error budget in [3] from 
the two-loop corrections and the higher-order nuclear-size effects. A conclu- 
sion of the present investigation is that it was originally underestimated [5] 
and actually it is approximately four times larger and thus gives the dominant 
effect to the uncertainty. 

The one-loop self-energy term can be presented in the form 

DlfiZ) = -iZarFsEiZ)EF, (5) 



where 

Fse{Z) 



021 In — + a2o + {Za)CsE{Z) 



(6) 



The coefficients a2i and 021 are known [6,7] (see, e.g., [3] for detail) 



021 = 7 ln2 = 3.303 22... , 
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020 = -5.221 23... (7) 
and evaluation of Cse{Z) is a purpose of our study. 

A perturbative evaluation of this coefficient has not yet been done. The results 
of [5] and this paper were obtained by fitting the numerical data obtained in [5] 
for a medium Z region, while the value of interest is related to Z — 1,2 . The 
calculation itself is a complicated problem, in fact much more complicated than 
the extrapolation. However, it is of a minor importance if the extrapolation is 
not performed properly. 

Some time ago we suggested an appropriate approach described in part in [8] . 
In fact, the authors of [5] stated that they applied our procedure and reached 



Cse(1)= 2.07 ±0.25 , 

Cse(2) =2.01 ± 0.19 . (8) 

The idea of the procedure, which we present here in more detail, is based on 
acknowledging a few typical features of any extrapolation for the self-energy 
contribution. 
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Fig. 1. Parabolic extrapolation of numerical data taken from [5] and its perturbation 
by various logarithmic functions. The open squares are the results of the original 
extrapolation in [5]. 

• The series contain numerous logarithmic contributions. 

• The data are available for medium Z (say from 10 to 30). 

• The slowly-changing logarithmic part (ln(Z)) can be hardly separated from 
the nonlogarithmic terms, which may include ln(Q;). 

• Because of that the strategy of the fitting it to exclude any logarithmic 
dependence on the first stage and perform a polynomial fit, e.g. C(Z^ = 
Aq + AiZ + A^Z"^ . Since each Z is accompanied with a, the coefficients 
decrease with their order S> 3> /I2. 

• The uncertainty of the fit is determined by stability of the results against a 
perturbation by logarithmic terms, e.g. by a difference with the results from 
fits C{Z) = Ao + A\Z + A2Z'^ + ACiog, where ACiog are various possible 
logarithmic contributions with coefficient of the natural value. 

However, a crucial part of the procedure is estimation of values of coefficients 
for the logarithmic terms which is explained in detail below while applying to 
the one-loop self energy contribution to D2\. 

The most general expression, which includes all contributions to Cse(-Z^) up 
to order {ZoLf is of the form 



The inclusion of higher order terms is not necessary here, since it will not 
improve accuracy of the fitting over Z < 30. 

The central value of the extrapolation according to our procedure can be 




(9) 
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obtained by setting 042 = 041 = 051 = 0. That is related to a parabolic fit (see 
Fig. 1) and on this issue we agree with results from [5]. However, we strongly 
disagree on the uncertainty of the result, i.e. we estimate a possible shift of 
such a result in a different way by using the logarithmic perturbation. 

A crucial question while estimating such a shift properly is related to the 
natural value of the coefficients. It is worth noting that odd and even terms 
of the series have different structure and magnitude of the coefficients. With 
this observation we estimate 

042 =±^21 — ±3.30 , 



where we apply a preliminary value of 030 from the parabolic extrapolation, 
which is sufficient for this purpose. Here a factor of vr appears because that 
is a characteristic value of the constant beyond the logarithm. The natural 
values of the unknown coefficients are pretty large and their underestimation 
should lead to a serious overestimation of the accuracy of the extrapolation 



While the introduced higher order corrections are negligible for Z — 1,2, they 
are not small for Z — 10 — 30, i.e. in the range of the input data for the 
extrapolation. The data [5] for Z below 10 look correlated and uncertain and 
we do not use them for the fitting. The same was done in [5]. 

Results of our extrapolation are summarized in Fig. 1 and Table 1. A scatter 
of the extrapolation results around the parabolic values allows us to estimate 
the uncertainty and our final results 



are approximately four times less accurate than in (8). Since we consider log- 
arithmic coefficients of order of unity in the natural units, any more accurate 
result may be achieved only on base of additional information and we consider 
the uncertainty in [5] as unappropriate. 

We have to comment briefly the stabihty of the fits. The parabohc coefficients 
a42 and 041 look unstable, however, one has to understand, that they are not 
related to any 'true' coefficient in (9). They arc 'effective' coefficients. The 
logarithmic terms ln(Z) slowly depend on Z and effectively the effective non- 
logarithmic coefficients include In(Zo), where Zq ~ 20 is kind of a medium Z 
value. If we like to see stability of the non-logarithmic coefficients, we have to 




(10) 



C{Z) to Z= 1,2 (cf. [5]). 



Cse(1)=2.1±1.0, 
Cse(2)=2.0±0.7 



(11) 
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Tabic 1 

Parameters of the parabolic extrapolation (0) and its perturbation by the logarith- 
mic functions (1-6). The logarithmic parameters are introduced as a fixed pertur- 
bation, while the non-logarithmic terms are found via least-square minimization. 
The uncertainty shown is a pure statistical one. 
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a42 





3.3 


-3.3 














041 











10.38 


-10.38 








040 


6.1(6) 


-1.98(6) 


-10.2(9) 


12.5(6) 


-24.7(8) 


-6.8(6) 


-5.4(7) 


a5i 

















2 


-2 


^50 


2.4(1.8) 


-13.1(1.7) 


17.9(2.4) 


-28.4(1.7) 


33.2(2.3) 


1.8(1.7) 


3.0(7) 


Cse(1) 


2.05(5) 


3.00(5) 


1.09(7) 


2.66(5) 


1.43(7) 


2.08(5) 


2.01(5) 


Cse(2) 


2.00(5) 


2.71(5) 


1.29(7) 


2.48(5) 


1.53(6) 


2.03(5) 


1.97(5) 



perturb the parabolic fit by ln(Z) — In(Zo). Such a fit would be more stable, 
however, mathematically the procedure is equivalent to what we have done. 
The actual stability should be seen not through values of auxiliary fitting 
parameters, but as a distribution of the results of the extrapolation Cse(1) 
and Cse(2). 



Wc have calculated above the one-loop self-energy contribution for hydrogen 
and deuterium {Z = 1) and heluim-ion {Z = 2). The latter is of most interest 
since it is a more sensitive test of bound state QED (see [3] for detail). 

Finally, we obtain for the helium ion 

Df^''\^B.e+) = -1 190.08(15) kHz , (12) 

while the most accurate experimental result is [9,10] 

(^He+) = -1 189.979(71) kHz . (13) 

A comparison of theory and experiment is summarized in Fig. 2, where the 
experimental data are labelled with the date of measurements of the hyperfine 
interval of the 2s state: 1958 [11] and 1977 [10]. Meanwhile the Is hyperfine 
interval was obtained experimentally in [9]. 

Theory is in perfect agreement with experiment. The most important conse- 
quence of our conservative estimation of the theoretical uncertainty is the fact 
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Fig. 2. The D21 difference for the hehum ion ^He"^: theory versus experiment. 

that the experimental result (13) is now twice more accurate than that of the 
theoretical prediction (12), while previously the relation was opposite. 

The procedure presented here has allowed to estimate uncertainty of the ex- 
trapolation of the numerical data properly for D21. Such an extrapolation 
occurs for a number of problems in QED theory for hydrogen-like atoms, and 
a plausible estimation of the accuracy is a crucial issue for the comparison of 
precision theoretical and experimental results. We hope that an application of 
our procedure will be helpful for other problems. 

The hydrogen and deuterium data will be discussed elsewhere together with 
experimental progress in the field. 
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